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ABSTRACT 

Detailed heat transfer measurements and predictions are 
given for a power generation turbine rotor with 129 deg of nom- 
inal turning and an axial chord of 137 mm. Data were obtained 
for a set of four exit Reynolds numbers comprised of the de- 
sign point of 628,000, -20%, +20%, and +40%. Three ideal 
exit pressure ratios were examined including the design point of 
1.378, —10%, and +10%. Inlet incidence angles of 0 deg and 
±2 deg were also examined. Measurements were made in a lin- 
ear cascade with highly three-dimensional blade passage flows 
that resulted from the high flow turning and thick inlet bound- 
ary layers. Inlet turbulence was generated with a blown square 
bar grid. The purpose of the work is the extension of three- 
dimensional predictive modeling capability for airfoil external 
heat transfer to engine specific conditions including blade shape, 
Reynolds numbers, and Mach numbers. Data were obtained by a 
steady-state technique using a thin-foil heater wrapped around 
a low thermal conductivity blade. Surface temperatures were 
measured using calibrated liquid crystals. The results show the 
effects of strong secondary vortical flows, laminar-to-turbulent 
transition, and also show good detail in the stagnation region. 


NOMENCLATURE 

A heater area [m 2 ] 

C x blade axial chord [mm] 
d leading edge diameter [mm] 

h heat transfer coefficient [W/m 2 *K] 

Fr Frossling number. Ft = Nu • ( d/C x )/ y/RelJn 
I heater current [Amps] 

K pressure gradient parameter, K — (pf pU 2 )(dU/ds) 
k thermal conductivity [W/m-K] 


k+ equivalent roughness height 
M Mach number 
Nu Nusselt number 
P pressure [Pa] 

Pr Prandtl number 

PR exit pressure ratio, PR = P' n /P €X 

q n heat flux [W/m 2 ] 

r recovery factor, r — Pr 1 / 3 

Re ex Reynolds number, Rc ex = pU ex C x /p 

s blade surface coordinate, normalized by blade span 

T temperature [K] 

T u turbulence intensity 

U total velocity [m/s] 

V heater voltage [Volts] 

t/ + equivalent normal distance 

z spanwise coordinate, normalized by blade span 

a incidence angle [degrees] 

7 specific heat ratio, 7 = 1.4 

e emissivity 

A x longitudinal integral turbulence length scale [mm] 

// dynamic viscosity [kg/s-m] 

a Stefan-Boltzmann constant 

Subscripts 

aw adiabatic wall temperature 

i c liquid crystal 

€X exit isentropic value 

i n inlet freestream value 

i s isentropic value 

Superscripts 
f total conditions 
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INTRODUCTION 

Turbine blade aerodynamics is generally given a high 
degree of analysis prior to commitment for fabrication. This 
attention to detail is carried through to the manufacturing 
specifications and quality inspections in the form of dimen- 
sional accuracies demanded of both the airfoils and the rotor 
disks. While a substantial amount of data exists concern- 
ing the effects of blade incidence angle, Reynolds number, 
and pressure ratio (loading) on aerodynamics, less informa- 
tion is available on the associated airfoil heat transfer dis- 
tributions for off-design conditions. Computational Fluid 
Dynamics (CFD) is increasingly being relied upon in the 
design and analysis of gas turbine components. The need 
to predict heat transfer along with aerodynamics during 
the design of turbine blades greatly complicates these anal- 
yses. Thus, the heat transfer predictive capability of CFD 
currently lags that of aerodynamics. Heat transfer data is 
therefore needed both to assess the effects of various flow 
parameters and to improve CFD analyses so that these ef- 
fects can be accurately predicted. 

Arts et ah, (1998) studied incidence, Reynolds num- 
ber, and exit Mach number effects on heat transfer, but 
with midspan thermocouples only. Incidence was varied 
over a range of —14 deg to +11 deg. Camci and Arts 
(1991) also examined incidence effects on a film-cooled heat 
transfer test blade, but again only at midspan. The three- 
dimensional nature of the flow and heat, transfer indicates 
that a full-blade understanding is needed, not just a two- 
dimensional understanding at midspan. Previous studies 
have shown that rotor geometries in linear cascades provide 
good midspari data as compared to their rotating equiva- 
lents. See, for example, Baughn et al. (1995) or Guenette 
et ah (1989). Also, Graziani et ah (1980) studied the effects 
of thick and thin endwall boundary layers on rotor heat 
transfer in a low speed linear cascade. Blair (1994) mea- 
sured heat transfer in a low speed facility on a rotating blade 
that had the same midspan section as the blade of Graziani 
et ah (1980). Similar three-dimensional heat transfer pat- 
terns were observed at and below midspan of those two 
tests, validating the use of linear cascades for these types 
of heat transfer studies. Dunn et ah (1994) measured rotor 
heat transfer at discrete chord wise and spanwise locations 
in a shock tube facility at engine-typical gas-to-wall tem- 
perature ratios and Mach numbers, These measurements 
were on engine-size hardware using heat flux gauges. Con- 
sequently, the resolution was less than for the large scale 
facility measurements. 

Detailed heat transfer data at conditions approximating 
those in actual engines are needed to assess flow parameter 
effects and to verify CFD predictive analyses. Variations in 
flow parameters often have an impact on secondary flows 
and thus on the three-dimensional nature of the heat trans- 


Parameter 

Value 

axial chord 

137.4 mm (5.410 in.) 

pitch 

130.0 mm (5.1 19 in.) 

span 

152.4 mm (6.000 in.) 

d, leading edge 

10.6 mm (0.417 in.) 

turbulence grid 

25.4 mm square bar 

blade passages 
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(d = 2 x minimum radius of cu nature) 



Figure 1. Test Blade Geometry 


fer distribution. Data is therefore needed all over the blade 
surface, not just at midspan. 

The blade geometry of the present study is represen- 
tative of an aerodynamically upgraded first stage turbine 
blade inidspan section for an established G.E. power gen- 
eration turbine. In an effort to continuously improve tur- 
bine performance, output, and durability, older designs may 
be upgraded in a number of ways, including aerodynam- 
ics, cooling, materials, and sealing. The present design is 
used in a non film-cooled application for which knowledge of 
the external convective heat transfer distribution is desired, 
both for design validation and improvement of predictive 
capabilities. The turbine blade is for a 2000° F class ma- 
chine. The midspan section design operating condition for 
full power yields an airfoil exit Reynolds number of 028,000 
with an inlet Mach number of 0.41 and an inlet total to exit 
static pressure ratio of 1.378, corresponding to an isentropic 
exit Mach number of 0.093. The airfoil section is entirely 
subsonic and has a total turning of 129 deg with a design 
inlet flow angle of 61.3 deg. The blade geometry is shown 
in Fig. 1 along with geometric details. Surface distance co- 
ordinates are also shown on the blade. 

The present study obtains both airfoil pressure and 
heat transfer coefficient distributions for a parameter design 
space which represents a set of limits that could be imposed 
by design uprate targets, manufacturing cost reductions, or 
operational demands. Incidence angle has been altered from 
-2 to +2 deg, exit Reynolds number from ”20% to +40%, 
and pressure ratio from -10% to +10%. These pressure ra- 
tio variations correspond to isentropic exit Mach numbers 
ranging from 0.563 to 0.794. The set of data produced in 
this study serves as a test of CFD predictive capability for 
an actual airfoil design space. The primary objective of this 
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Figure 2. Transonic Turbine Blade Cascade Facility 


research is the continued extension of three-dimensional pre- 
dictive modeling techniques for airfoil external heat transfer 
through the acquisition of full-surface blade data. The air- 
foil and flow conditions of the present study model those 
of an engine design, including a moderate level of surface 
roughness as provided by the test method. The airfoil shape 
is distinctly different than one previously tested in the same 
facility, thus providing a new CFD test case for previously 
validated predictive techniques. The CFD comparisons are 
challenged further by the additional variation of incidence 
angle and by the relatively low Reynolds numbers. The 
CFD predictions, when compared to the measurements, 
point out areas where modeling improvements are needed. 
Liquid crystal surface temperature measurement techniques 
are capable of obtaining detailed data over an entire surface 
and were therefore chosen for the present study. 


DESCRIPTION OF FACILITY 
Cascade Description 

The NASA Glenn Research Center Transonic Turbine 
Blade Cascade provides heat transfer and aerodynamic data 
to verify CFD analyses. The facility can and has been used 
to test over wide, independent ranges of Reynolds number 
and Mach number. High pressure, ambient temperature 
air can be throttled to an inlet pressure that can be varied 
from the exit pressure to over one atmosphere. The facility 
is designed for independent control over exhaust pressure 
which can be varied from the inlet pressure down to 15.9 kPa 
(2.3 psia). A two-dimensional view of the facility is shown 
in Fig. 2. The heat transfer measurement blade was the 
sixth blade from the left, and the pressure measurement 
blades were the fifth and sixth. 

While the blade and cascade geometries are two- 
dimensional, the blade passage flows are highly three- 
dimensional due to the turning and to the thick boundary 
layers developed on the long cascade inlet. The thickness 


of each endwall boundary layer was approximately 35% of 
the half-span, as measured approximately one axial chord 
upstream of the blade row. The thick boundary layers and 
high flow turning resulted in spanwise flow variations as 
large as those seen in rotating turbines. Sec, for example, 
Joslyn and Dring (1992) or Tliulin et al. (1982). These vari- 
ations result from secondary flows, and accurately predict- 
ing them is a significant test, of a three-dimensional analysis. 

Inlet turbulence was produced using a square bar, bidi- 
rectional grid upstream of the blade row. The grid was 
made of 25.4 mm (1 inch) square hollow bars with 75 uni- 
formly spaced 3.2 mm (1/8 in.) diameter holes blowing air 
in the upstream direction. One bar extended between the 
inlet boards at midspan and three bars spaced six bar di- 
mensions apart were placed spanwise as shown in Fig. 2. 
The total mass flow from the grid was 0.68 kg/s (1.5 lb m /s) 
for all cases, corresponding to 11% of the cascade mass flow 
at the design exit Reynolds number. The open area of the 
grid was 74% and it was located 1.10 m (43 in.) upstream of 
the heat transfer measurement blade. It was positioned nor- 
mal to the inlet flow as shown in Fig. 2. An aerodynamic 
probe measurement slot was located an axial distance of 
127 mm in front of the cascade face. The inlet turbulence 
level measured there was approximately 9%. Using the cor- 
relation of Baines and Peterson (1951), this can be expected 
to decay to about 7.5% at the cascade face. The longitu- 
dinal, integral turbulence length scale was measured to be 
A a- = 22 mm. Pitchwise and spanwise surveys of mean ve- 
locity, Tu , and length scale were made in the probe slot over 
several blade passages. The results of those measurements 
can be found in a previous study by Giel et al. (1999). The 
blades of the current study have the same pitch and are 
located at the same pitchwise locations as those of the pre- 
vious study. 

Figure 2 shows the test section mounted on a large ro- 
tatable disk. The disk can be rotated H-15 deg to -30 deg 
to accommodate a blade with a different inlet flow angle 
than the original as well as to vary the incidence angle. Up- 
stream inlet boards were used, but the facility uses no exit 
tailboards. Endwall static pressure measurements and other 
flow field aerodynamic probe measurements all showed ex- 
cellent periodicity in at least the center three blade passages 
(see Giel et al., 1996). 


Pressure Measurement Blade Description 

The two blades forming passage 5 were instrumented for 
pressure measurements. Blade 5 had 20 static pressure taps 
at 10% span and 20 more at 25% span. Blade 6 had 40 taps 
at midspan. The tap diameters were 0.5 mm (0.020 in.). 
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MEASUREMENT TECHNIQUES 


Heat Transfer Measurement Blade Description 


The blade in position 6 was fabricated from low conduc- 
tivity foam for heat transfer measurements. The undersized 
foam core was placed in a mold and gelcoat epoxy was cast 
around it. The outer epoxy shape was itself undersized by 
152 //m (0.006 in.) to accommodate a 25 fim (0.001 in.) In- 
conel thin foil heater and a 127 fim (0.005 in.) double-faced 
adhesive layer. The single sheet heater covered the entire 
blade except on the trailing edge circle. Two thin-film ther- 
mocouples were glued to the heater sheet at z — 0.25, one 
on the pressure surface at s — -0.29 and one on the suction 
surface at s = 1.51. These were used for overheat control 
and for liquid crystal temperature verification. Flat black 
paint was sprayed over the outer surface and a 25 x 25 mm 
(lxl in.) grid of white dots was painted on for location 
reference. Finally, chiral nematic, micro-encapsulated liq- 
uid crystals were sprayed onto the outer heater surface. A 
temperature calibration plate, instrumented with a thermo- 
couple, and a roughness measurement plate were sprayed at 
the same time. The liquid crystal yellow-line temperature 
was calibrated with the calibration plate and verified dur- 
ing data acquisition with the two blade surface-mounted 
thermocouples. Both ends of the measurement blade were 
made of 12.7 mm (0.50 in.) thick nickel-plated copper, cut 
out in the shape of the blade, and were used as buss bars 
to supply DC electrical power to the heater. The heater 
foil extended over the copper ends and was tack welded 
to them. Hollow studs extended from the copper, through 
the acrylic windows, and were used for mounting and for 
electrical connections. They were also connected through 
vacuum tubing to the exhaust section of the facility to ap- 
proximately equalize the pressure inside and outside of the 
blade, preventing damage to the surface when the test sec- 
tion pressure was changed. Both end walls were G3.5 mm 
(2.50 in.) thick clear acrylic for optical access. An outline 
of the window is shown in Fig. 2. 


A profilometer was used to measure several 5-mm traces 
on the roughness measurement plate after the liquid crystals 
were applied. A 60 x photomicrograph showed that the pro- 
filometer stylus did not alter the crystals. The arithmetic 
mean roughness was 6.5 f,im and the r.m.s. roughness was 
7.8 fin l. The average peak-to-peak distance was 97 fim. The 
equivalent roughness height was estimated to be 7 fim. The 
maximum normalized roughness height, fc + , was estimated 
to be 1.2, so from a surface roughness standpoint the blade 
surface was hydraulically smooth. This tested condition is 
thus not typical of in-service operating blade roughness. 


Blade Static Pressure Measurements 

Surface static pressures were measured independently 
from the heat transfer measurements using a separate blade 
set. The two blades forming passage 5 were instrumented 
with 80 static pressure taps. The pressure taps were 
scanned at a rate of once per second with an electronically 
scanned measurement system. Fifteen scans were averaged. 
The inlet total pressure was measured as the average of 
three midspan Kiel probes upstream of passages 4, 7, and 
8. The probes were located an axial distance of 127 mm 
(5.0 in.) ahead of the cascade face. More details of the 
blade static pressure measurement techniques are given by 
Giel et al. (1996). 

Heat Transfer Measurements 

The selections of liquid crystal yellow line tempera- 
tures and color band widths were made subject to several 
constraints. The crystal temperature was chosen to be as 
high as possible to minimize uncertainty without exceeding 
the 80°C (175°F) material limit of the double-faced adhe- 
sive. This limit was approached in regions of minimum heat 
transfer when the heater power was increased to make mea- 
surements in regions of maximum heat transfer. A crystal 
temperature 14° C (2o°F) above the inlet total temperature 
typically met both criteria. Because the inlet temperature 
varied from day to day, two crystals were mixed and sprayed 
onto the surface. Data were then obtained from whichever 
crystal offered a higher AT without exceeding the material 
limits. Often, a combination of crystal temperatures were 
used to map the entire surface. Narrower band width cr} r s- 
tals provide better resolution in regions of very low gradi- 
ents but the yellow line could become too thin to be clearly 
visible in high gradient regions. Crystals with full color 
bandwidths of ±2°C were found to be a good compromise 
for the current measurements. 

The following procedure was used to obtain the heat 
transfer data: Flow conditions were established and the 
heater power was increased until the first yellow isotherm 
was visible. When steady state conditions were achieved, 
typically within 5 minutes, data recording of voltage, cur- 
rent, and flow conditions was initiated. The blade surface 
was simultaneously photographed with four 35 mm still 
cameras using color slide film . Each camera photographed 
a subset of the blade surface, with some overlap between 
views. High speed photographic strobes were used to pre- 
vent radiative heating of the liquid crystals. The strobes 
were mounted as close as possible to the camera lens to 
minimize differences between viewing and illumination an- 
gles. Moffat (1990) showed that micro-encapsulated crys- 
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tals are less sensitive to illuminating and viewing angles 
than non-encapsulated crystals. No differences in isotherm 
location were detected between overlapping photographs. 
Heater power was then increased to move the isotherms and 
the process was repeated until the entire blade surface was 
mapped. Eighteen power levels were typically used for each 
flow condition. 

The difference between the surface temperature and the 
local adiabatic wall temperature was used to define the heat 
transfer coefficient. In regions where data were available 
from both crystals, agreement between the crystals verified 
that the correct definition of h was used and that the adia- 
batic wall temperature was reasonably accurate. 

The following procedure was used to reduce the heat 
transfer data. The photographic slide images were projected 
onto paper and the reference dot pattern and the isotherms 
were sketched. The slide image dot pattern was aligned for 
subsequent isotherms. Separate sketches were made for each 
camera and crystal temperature. These isotherm sketches 
were then digitized, with between 3000 to 5000 digitized 
points for each flow condition. A grid having five times 
the resolution of the blade dot pattern painted was laid 
over the blades and photographed. Every fifth point on 
the grid photograph was aligned with the dots in the data 
photographs. The finer resolution coordinates were then 
used to interpolate each digitized point to (s, z) unwrapped 
blade coordinates. This procedure corrected for distortions 
due to blade curvature and viewing angles. 

The local surface heat flux, </', was determined from 
the heater voltage and current., then corrected for radiative 
heat loss: 


<7" = El - ea(TH - T*J 

With e = 0.98, radiative losses were at most 6% of the 
net heat flux and typically much less. The heat transfer 
coefficient and the Nusselt number were defined as follows: 


(Tic - T aw ) 


and 


Nu = 


h • C x 
k(T[ n ) 


The local adiabatic wall temperature, T au ,, is: 


(i) 


T na _ 1 ~ r 

TL 1 1 + 0.5(7 - 1 )Ml 

The choice of T aw as the convective driving temperature 
ensures that h and Nit are reasonably independent of the 


thermal boundary condition, specifically, independent of the 
particular choice of liquid crystal temperature. The inlet to- 
tal temperature, T/ rl , was determined by a mass weighted 
average of the main flow and the turbulence grid flow. Two 
additional total temperature probes were located an axial 
distance of 127 mm (5.0 in.) ahead of the cascade face 
and agreed with the mass weighted average to within the 
measurement uncertainty. The local adiabatic wall temper- 
ature was needed at each digitized point in order to cal- 
culate the heat transfer coefficient. This information was 
interpolated from CFD calculations of blade surface static 
pressure which will be shown to agree well with measured 
values. The calculated values provided significantly better 
spatial resolution than was available from the experimental 
data. The isentropic Mach number, Mj Si was determined 
from CFD calculations, and a recovery factor of r = Pr 1 / 3 
was used everywhere. The data were triangulated for con- 
tour plotting and interpolated onto constant span lines of 
15%, 25%;. and 50% for line plotting. 


Uncertainty Analysis 

An uncertainty analysis was performed using the 
method of Kline and McClintock (1953). The major 
sources of uncertainty are the adiabatic wall temperature. 
6T aw — ±1.0°C (=tl.8°F) and the liquid crystal tempera- 
ture, 8Tic = d=0.5°C (±0.9°F). At the temperature differ- 
ences and the relatively low Mach numbers of the current 
study, the uncertainty in pressure measurements, the differ- 
ences between measured and calculated pressures, and the 
interpolations are minor contributors to the overall uncer- 
tainty of _ Taw as compared to the uncertainty in T' n . The 
uncertainty in T aw for the same reason assumes no uncer- 
tainty in the recovery factor, r, although it is not precisely 
known as shown by Schlichting (1979). Other less domi- 
nant sources of uncertainty such as the heater voltage and 
current, SV/V = 81 /I = 1 %, and the heater length and 
width, 5L/L — 2% and 8W/W = 1 %., were included in the 
analysis. The overall uncertainty in Nu was determined 
to be less than 13%; for all cases and was typically much 
less. Because the uncertainty varied with specific condi- 
tions, maximum values for each case will be presented with 
the results. 

Another possible source of bias uncertainty is ther- 
mal conduction within the measurement blade. A two- 
dimensional finite-difference conduction analysis w'as per- 
formed to examine this error. A typical, measured midspan 
h distribution was prescribed on the outer surface, and tem- 
peratures within the blade were calculated. A thermal con- 
ductivity of 0.2 W/ni • K was used for both the foam and 
the gelcoat epoxy. Conduction within the blade tended to 
increase q" in cooler regions and decrease it in hotter re- 
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gions. Corrected values of h were then determined from 
Eq. (1). Figure 3 shows that the maximum conduction er- 
ror is approximately 6%, which is within the overall uncer- 
tainty. To assess the validity of the conduction analysis, the 
finite- difference grid resolution was doubled in the leading 
edge region. The finer grid results were essentially identical 
to the results shown in Fig. 3. 


CFD ANALYSIS 

To illustrate where analysis improvements are needed, 
baseline CFD results were done using the three-dimensional 
Navier-Stokes analysis code, RVC3D. described by Chima 
and Yokota (1990) and by Chima (1991). C-type grids were 
generated using the method of Arnone et al. (1992). Pre- 
dictions were made using a 313 x 49 x 49 grid with a near 
wall spacing, y+ < 1. 190 grid lines were wrapped around 
the blade. A uniform wall temperature of 1.07 x T( n was 
specified on the heated blade and a uniform wall temper- 
ature of 1.0 x T’ n was specified on the unheated end wall. 
The flow over much of the blade was expected to be turbu- 
lent so little difference would be expected from uniform heat 
flux and uniform wall temperature heat transfer solutions 
(Kays and Crawford, 1980). This was verified by running a 
calculation with a wall temperature boundary distribution 
consistent with one that would be present for a uniform 
heat flux. The two solutions agree very closely. All of the 
analyses shown herein assumed midspan symmetry. The so- 
lutions were mirrored across midspan for easier comparison 
with the full-span experimental data. 

A two-layer algebraic turbulence model, described by 
Chima et al. (1993), was used. The model was designed 
to avoid spurious inner-to-outer region crossovers, which 
sometimes occurs with Baldwin-Lomax type algebraic tur- 
bulence models. The inner and outer blending modification 
described in Giel et al. (1999) was again used. An algebraic 
model, Smith and Kuethe (1966), to account for the ef- 
fects of freestream turbulence on laminar heat transfer was 
used. This model was found by Giel et al. (1999), and by 
Boyle and Simon (1999) to improve the agreement of pre- 


dicted heat transfer with the experimental data. However, 
for the results presented herein this model did not improve 
the agreement with data. Predictions are shown with and 
without this laminar augmentation model to illustrate the 
significance of the model for heat transfer predictions. 

Transition start was specified using Mayle’s (1991) 
model. Tu at transition start was the upstream Tu times 
the ratio of the upstream velocity to the local isentropic ve- 
locity. This approach was found by Boyle and Simon (1999) 
to give an accurate prediction for the start of rotor blade 
transition. The transition length model of Boyle and Si- 
mon (1999) was used. This model is an extension of the 
one by Solomon et al. (1995) to include Mach number ef- 
fects. 


BLADE STATIC PRESSURE RESULTS 

Figure 4 shows sample measured and calculated isen- 
tropic surface Mach number distributions on an unwrapped 
blade surface. Fig. 4(a) shows calculated contours of Mj s 
for the design case. Figs. 4(b-d) show comparisons with 
data for variations in PR. In all cases, the comparison 
with data is excellent. The surface contour plot shows the 
strong three-dimensionality of the flow, as evidenced by 
the distorted suction surface contours caused by the horse- 
shoe/passage vortices. These vortices are seen to exit the 
blade row at about 25% span. The line plots also show that 
the 25%: span Mach number distribution exceeds that of the 
10% and 50% on the back half of the suction surface. No 
decelerating flow regions are seen on the suction surface un- 
til after the throat at s & 1.0 where very slight deceleration 
occurs on the uncovered portion of the blade. Deceleration 
is also seen on the pressure surface just downstream of the 
leading edge, extending to s & —0.4. These regions are 
noted because of their potential implications on transition 
and thus on heat transfer. The calculated values were used 
for heat transfer data reduction because they provided sig- 
nificantly better spatial resolution than was available from 
the experimental data. 


EXPERIMENTAL HEAT TRANSFER RESULTS 

Data and CFD calculations will be presented as Nus- 
selt number contours on the unwrapped blade surface and 
as line plots at 15%, 25%, and 50% span. Discussion on the 
experimental data all will be presented first, followed by a 
discussion of the calculated results. Gray areas on the data 
contour plots indicate regions where data were not available. 
The copper buss bars will also be evident in the data plots 
as regions near the endwalls where no data were available. 
The CFD calculations were made with a uniform temper- 
ature specified on the entire blade and did not model the 
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Figure 4. Measured and Calculated Blade Loadings 


buss bar thermal boundary condition. Lines showing their 
locations are included in the contour plots for reference. It 
should be noted that no midspan symmetry was assumed 
in presenting the experimental data; measured symmetry is 
typically excellent, but minor differences can be seen. The 
first case to be presented will be the baseline case, obtained 
at the design flow conditions. Subsequent cases are grouped 
to examine the effects of Reynolds number, pressure ratio, 
and incidence angle. 

The measured flow conditions for the eight heat transfer 
cases are listed in Table 1. Repeatability of Re ex was within 
1% for all cases and repeatability of PR was within 0.2%. 
As mentioned earlier, Nu uncertainty depended primarily 
on the difference between the air temperature and the liquid 
crystal temperature. The variations in uncertainty shown 
in Table 1 resulted from day-to-day variations in inlet air 
temperature. Red, in is included in the table for discussion 


Table 1. Description of Blade Heat Transfer Cases 


Case 

Re rx 

PR 

a in 

Redan 


peak 

Fr 

max Nu 
uncert 

i 

G22 : 100 

1.376 

0.0° 

29,939 

0.371 

0.828 

9.8% 

2 

500.800 

1.381 

o 

p 

24,271 

0.376 

0.950 

7.6% 

3 

751. GOO 

1.37G 

0.0° 

35,945 

0.369 

0.866 

7.5% 

4 

880.400 

1.378 

0.0° 

42,151 

0.370 

0,869 

10.8% 

5 

626,900 

1.240 

0 

O 

o 

30,801 

0.327 

0.839 

12.4% 

6 

625,700 

1.513 

0.0° 

29.068 

0.389 

0.998 

7.7% 

7 

G2G.700 

1.377 

-2.0° 

27,498 

0.335 

0.899 

9.9% 

8 

629.500 

1.376 

+2.0 C 

33,250 

0.414 

p 

oo 

•q 

7.4% 


of the stagnation point heat transfer. An inlet Reynolds 
number based on axial chord can be calculated from this by 
multiplying by C x jd— 13.0. 

Baseline Case 

Figure 5 shows the measured and predicted Nusselt 
number distribution at the design flow condition of Re ex = 
628,000, PR = 1.378, and a = 0 deg. The pressure side 
data shows a very flat, two-dimensional Nu distribution 
with a slight increase at the trailing edge, possibly indicat- 
ing the start of transition. A leveling of Nu near s = -0.1 
can be seen in Fig. 5(c). Just downstream of the stagnation 
point, the data shows no indication of transition, despite 
the high freestream turbulence intensity and the decelerat- 
ing flow region, —0.4 < s < —0.1. The heat transfer data, 
however, may not be a good indicator of transition. The low 
velocities could counteract the effects of transition, leading 
to uniform Nu levels. This will be verified with comparisons 
at different Reynolds numbers. 

The stagnation region is clearly evident in the figure, 
with a somewhat steeper gradient on the pressure side than 
on the suction side. Table 1 lists peak Frossling numbers, 
defined as Fr = Nu • ( d/C x )/ \/ Red, in* The peak Nusselt 
numbers were corrected for blade conduction by a factor of 
1.06 as suggested by Fig. 3. This correction was only ap- 
plied to the peak data of Table 1 where the conduction error 
was largest, not to the the remainder of the data where it 
has less effect. An ellipse with an aspect ratio of 1.5 was 
found to best describe the blade leading edge. Van Fos- 
sen et al. (1995) studied stagnation region heat transfer on 
isolated elliptical leading edges. They indicate that a lami- 
nar Frossling number of 0.870 is appropriate for this aspect 
ratio. They also developed the following correlation for a 
stagnation point turbulent augmentation factor: 


F / / 4 \ - 0.574 

= 0.0085 \ Tu Re° d s ( ) +1.0 (2) 

Friam V V d ) 
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Figure 6. Heat Transfer Distributions - Case 2 

Despite the relatively high Tu of the current study, this cor- 
relation predicts an augmentation factor of only 1.12, and 
thus a Frossling number of 0.976. The low augmentation 
results from the low Reynolds number and high A x /d ra- 
tio. The measured value of 0.828 is about 15% less than 
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Figure 8. Heat T ransfer Distributions - Case 4 

the correlation predicts. This could be due to an under- 
prediction of conduction effects or fundamental differences 
between heat transfer rates on isolated ellipses at zero angle 
of attack and a cascade of highly loaded blades. 

Just downstream of the stagnation region on the suc- 
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(a) measured Nu contours 





tion surface, a very slight increase in Nu is observed near 
0*2 < s < 0.4. This could indicate that transition started, 
but was damped by the accelerating flow. Farther back 
on the suction surface, the effects of the strong secondary 
flows are evident. Many of these effects can be attributed 



(a) measured Nu contours 




(a) measured Nu contours 




to the vortex structures described by Langston (1977). The 
passage vortex and the pressure-side leg of the horseshoe 
vortex approach the suction surface. Heat transfer is en- 
hanced by the relatively cool secondary flow fluid impinging 
on the endwall regions of the suction surface. The vortices 
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Table 2. Effects of Reynolds No. on Pressure Surface Heat Transfer 






.S’ = 

-0.26 

s = 

-0.51 

Case 

Re„ 

ANii( am 

A A Ufurb 

Nu 

ANu 

Nu 

ANu 

1 

628,000 


— 

712 


675 

- 

2 

—20% 

-11% 

-16% 

603 

-15% 

607 

-10% 

3 

+20% 

+ 10% 

+ 16% 

817 

+ 15% 

751 

+ 11% 

4 

+40% 

+ 18% 

+31% 

921 

+29% 

816 

+21% 



lift off the endwall and approach midspan, exiting the rotor 
at around 25% span. Their effect is seen in Fig. 5(c) as 
departures from midspan data; first at 15% span, s = 0.7; 
then at 25% span, s = 1.1. The data of Fig. 5(c) shows 
that midspan transition occurs over the approximate range 
0.7 < s < 1.1. Midspan Nu increases from G25 at s = 0.7 
to 824 at s — 1.4. Note that this transition start occurs 
before the slight decelerating flow region. At 25% span, 
secondary flows enhance the heat transfer for s > 1.0. This 
enhancement makes a determination of the end of transition 
difficult. 

Effects of Reynolds Number 

Table 2 lists and shows comparisons of Nusselt num- 
ber data at two midspan points on the pressure surface. 

i /f> 

&Nuiam. was determined by assuming Nu oc ReFF and 

A t T 

ANuturb was determined by assuming Nu oc Re ex \ Note 
from Fig. 4 that the flow at s = —0.26 is decelerating while 
the flow at s = —0.51 is strongly accelerating, but that the 
velocities at —0.51 are quite low. 

Figure 6 shows results obtained at the design pressure 
ratio and inlet flow angle, but at a Reynolds number 20% 
less than design (Case 2). No transition start was observed 
near the end of the pressure surface as was seen for Case 1. 
The level region of Nu near s = -0.1 is still seen for this 
case. A comparison of heat transfer rates can be made at 
the midspan locations listed in Table 2. At s — —0.26 the 
change in Nu with Rc ex is consistent with expected changes 
in a turbulent boundary layer. At s = -0.51, however, 
the change in Nu with Re €X is consistent with expected 


changes in a laminar boundary layer. This indicates that 
the flow transitioned in the decelerating flow region and 
then relaminarized under strong acceleration. The pressure 
gradient parameter of I\ = 12 x 10 -6 at s ~ —0.51 supports 
this. K was greater than 3.5 x 10~ 6 for s < -0.42. 

The laminar stagnation point Frossling number of 0.870 
from Van Fossen et al. (1995) is independent of Reynolds 
number. Equation (2) predicts an enhancement of 1.11 and 
thus a Frossling number of 0.968. The measured value of 
0.950 agrees much better than for Case 1. 

Heat transfer on the suction surface appears to be quite 
similar to that of Case 1. except that the lower Reynolds 
number has delayed transition all the way to the throat 
and to the decelerating flow region. No increase in Nu is 
observed near 0.2 < $ < 0.4 indicating that transition never 
started in that region. Fig. 6(c) shows the midspan transi- 
tion region to be 1.00 < s < 1.35. The fact that transition 
occurred in a decelerating flow region would account for the 
doubled relative Nu rise in going from the laminar value of 
535 at s — 0.9 to the fully turbulent value of 864 at $ = 1.4. 
Quarter-span transition appears to start at the same loca- 
tion but also to be more abrupt and intense. At 15%: span, 
secondary flows apparently incited much earlier transition 
(. 8 as 0.6)'. 

Figure 7 shows the effects of an increase in Reynolds 
number 20%; over the design value. Pressure surface tran- 
sition is quite evident for 5 < -0.7. Table 2 shows that 
the changes in Nu with Re ex are again consistent with tur- 
bulent flow at s — —0.26 and with relaminarized flow at 
s — —0.51. The measured peak Ft — 0.866 is 12% lower 
than the correlation prediction, similar to Case 1. 

On the suction surface, definite increases in Nu are ob- 
served near 0.2 < s < 0.4 at 15%, 25%, and 50% span, indi- 
cating early transition despite the flow acceleration shown 
in Fig. 2. Nu levels gradually decrease after this, but then 
increase slightly to the throat at s ~ 1.0. At s = 0.7, 
Nu = 919, an increase of 47% over the Case 1 value at 
the same location. This is a much greater increase than 
would be expected if the flow had remained laminar to this 
point, again indicating that the flow of Case 3 was tran- 
sitional in this region. For $ > 1.0, midspan levels of Nu 
remain constant while quarter-span levels rise because of 
secondary flow effects. The 15% span results again show 
earlier transition and increased secondary flow heat trans- 
fer enhancement. 

Case 4 results are presented in Fig. 8. This is the high- 
est Reynolds number case, 40% over design. The trends 
noted in Case 3 are even more evident in this case. Pres- 
sure surface midspan transition starts at about the same 
location, s = —0.7. The flow appears to be fully turbulent 
near the end of the pressure surface. Table 2 again shows 
data consistent with turbulent flow at s — -0.26 and with 
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relaminarized flow at s = -0.51. At s = -1.0, the midspan 
Nusselt number is 1365 as opposed to 1168 for Case 3 at the 
same location. This 17% increase is in good agreement with 

I /c 

an expected increase of 13% assuming Nu oc Rc ex . This 
comparison could not be made for Cases 1 or 2 because 
the flow was either laminar or transitional at this location. 
The measured peak Fr — 0.869 is again 12%: lower than 
the correlation prediction, but agrees well with the previ- 
ous case. Transition is again evident near the beginning of 
the suction surface. The more pronounced decrease in Nu 
for 0.4 < s < 0.7 may indicate that complete transition had 
occurred by s ~ 0.4. At a midspan distance of s = 0.7, 
Nu increased from 919 for Case 3 to 1103 for Case 4. This 
20% increase is slightly more than the 13% expected in- 
crease from Reynolds number. This may further indicate 
that the flow in Case 3 was still transitional at this point. 
After the throat , midspan levels of Nu continue to decrease 
as the fully turbulent boundary layer grows, while quarter- 
span levels increase because of secondary flow effects. The 
15% span data is again significantly higher because of sec- 
ondary flow heat transfer enhancement. In comparing the 
15% span data from Cases 1-4 it is observed that whether 
or not the flow transitioned early as was the case for the 
higher Reynolds number cases, the data at 15%; span al- 
ways departed from the 25% and 50% data and increased 
sharply at s = 0.7. This gives an indication of the path of 
the horseshoe/passage secondary flow vortex. 

Effects of Pressure Ratio 

The isentropic exit pressure ratio was varied ±10% 
from the design value of 1.378. The independent inlet 
and exhaust control of the facility allowed the nominal exit 
Reynolds number to be held fixed at the design value of 
620,000 for these cases. Changes in the inlet Reynolds num- 
ber thus were minimal. Relative to the exit Reynolds num- 
ber, Re m changed only +2.3% for a -10%; change in PR 
and changed —3.3% for a -+10%: change in PR. 

The results for Case 5 are shown in Fig. 9. Compar- 
ing this to Case 1 in Fig. 5 shows the effects of an exit 
isentropic pressure ratio reduction of 10%. On the pressure 
surface, those effects are minimal. Fig. 9(c) still shows a 
slight increase in Nu at the end of the pressure surface, in- 
dicating the start of transition. Table 1 shows that the peak 
Frossling number is 0.839, again very close to that of the 
baseline case. Fig. 9(c) shows that the suction surface Nu 
distribution starts off the same and has the same qualita- 
tive trends, but ends up about 10% higher than the baseline 
case. No explanation for this can be found, except to note 
that because of day-to-day temperature variations, this case 
had the highest uncertainty, 12.4%. It is interesting to note 
that Arts et al. (1998) also showed nearly identical pressure 


surface behavior and gradually increasing heat transfer dif- 
ferences along the suction surface with decreasing pressure 
ratio. They attributed this to shock/boundary layer inter- 
action at their transonic Mach numbers, but it appeared to 
be evident even before the shock impingement location. 

Figure 10 shows the results for a 10% increase in PR. 
Comparing the contour plot of Fig. 10(a) and the 15%, 25%;, 
and 50% span line plots of Fig. 10(c) to those of the baseline 
case in Fig. 5 shows that the Nusselt number distributions 
are nearly identical. The only difference appears to be at the 
stagnation point where the Frossling number is 2% higher 
than the correlations would predict. This is in contrast to 
all of the other cases in which measured Frossling numbers 
were less that correlation predictions. No explanation can 
be offered for this. 

Effects of Incidence Angle 

Figure 11 shows the effects of —2 deg of incidence. The 
only appreciable difference from the baseline case appears to 
be on the pressure surface very near the stagnation region. 
Figure 5(c) showed a small, level region of Nu around s — 
—0.1. Figure 11(c) shows a slight dip in this same region 
that can be attributed to the negative incidence angle. This 
could indicate a small separation bubble or at least incipient 
separation. Arts et al. (1998) showed a similar, but much 
stronger effect of negative incidence for a -5 deg incidence 
case. 

Figure 12 shows the effects of +2 deg of incidence. The 
pressure surface and stagnation region are very similar to 
the baseline case, but the suction surface shows signs of 
early transition. The effects are similar to a 20% increase 
in Re ex as shown in Fig. 7. Note, however, from Table 1 that 
the inlet Reynolds number for Case 8 is 11% higher than 
the baseline case even though the exit Reynolds numbers are 
nominally the same. The early transition leads to increased 
Nu levels on the remainder of the suction surface. 


COMPUTATIONAL HEAT TRANSFER RESULTS 

Figure 5(b) presents computed Nu contours while 
Fig. 5(c) presents comparisons with data at 15%. 25%, and 
50% span. For this and all subsequent figures, calculation 
results shown in part (c) correspond directly to those shown 
in part (b), either with or without the laminar augmentation 
model. The analysis for this case does not agree well with 
the experimental data. The analysis methods used for the 
comparisons had previously been verified by comparisons 
with experimental data. Giel et al. (1999) showed compar- 
isons with experimental data using the same analysis, and 
the agreement was much better than for the results shown 
herein. The comparisons were for surface heat transfer on a 
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rotor geometry tested in the same facility and with the same 
turbulence grid as the present data. While the test Mach 
numbers, Reynolds numbers, and the rotor geometry all dif- 
fered from the current test, the blade row inlet conditions 
were similar. The blade of the current test has a signifi- 
cantly smaller leading edge diameter. The blade thickness 
just behind the leading edge is also significantly less for the 
current blade which results in lower mid-passage velocities. 
These both influence the scale and intensity of turbulence 
within the blade row. Modeling for transition and the effect 
of turbulence on pre-transition heat transfer were verified 
by comparisons with rotor and stator data as shown by 
Boyle and Simon (1999). They compared midspan predic- 
tions with data and showed good agreement . The modeling 
approach of the two-dimensional analysis was incorporated 
into the three-dimensional analysis. The areas of disagree- 
ment between the predictions and the data herein illuminate 
areas where improved modeling is needed. 

Figure 5 shows that the analysis significantly over- 
predicts the leading edge region heat transfer. The aug- 
mentation due to inlet turbulence is greater than measured. 
The model for augmenting leading edge heat transfer only 
accounts for turbulence intensity; it does not account for 
turbulence scale. As mentioned earlier, because of the large 
A x /d ratio, the expected leading edge augmentation was 
small. The leading edge Frossling number is close to 1.0 in 
the data, but the predicted value is nearly 1.5. The exper- 
imental data show transition beginning near s — 0.7. In- 
termittency information from the CFD code indicates that 
the predicted suction surface transition start occurred at 
about s = 0.25. Transition effects are not immediately 
evident in this region because the high heat transfer on 
the forward portion of the suction surface is due to aug- 
menting the laminar viscosity to account for the effect of 
freestream turbulence. This is the same model as was used 
in the leading edge region. Consequently, the heat trans- 
fer is over-predicted. The differences in the augmentation 
model predictions between the current study and that of 
Giel et al. (1999) is attributed to their larger circular lead- 
ing edge diameter-to-chord ratio and to their higher flow 
acceleration. 

On the uncovered portion of the suction surface, s > 1, 
the analysis under-predicts the midspan heat transfer, but 
over-predicts the heat transfer at 15% and 25% span. Be- 
tween 30 and 40% of span in this region, the predicted 
heat transfer changes rapidly. The temperature and/or the 
amount of fluid being swept up the suction surface are not 
adequately predicted. In this region, while the local pre- 
dictions are in error, the predicted span wise average heat 
transfer at each s location agrees reasonably well with the 
experimental spanwise average heat transfer. The analysis 
predicts the character of the heat transfer variation, but 


in the uncovered region, the variation is greater than the 
experimental variation. Additional calculations were per- 
formed with a heated end wall which showed, consistent with 
Giel et al. (1999), that the effect of the endwall thermal 
boundary condition on rotor heat transfer was relatively 
small, except very near the endwall. This indicates that 
the hulk of the secondary flow impacting the suction sur- 
face apparently came from the freestream. 

The pressure surface heat transfer predictions have 
common characteristics for all cases. Transition occurred 
in the forward portion of the pressure surface. Figure 4 
shows diffusion occurring between s = —0.1 and -0.4. The 
diffusion resulted in the analysis predicting transition. Nu 
decreases despite transition because the velocity decreases. 
For most of the pressure surface the predicted flow was 
turbulent. The pressure surface heat transfer predictions 
shown in Figs. 5 12 agree well with the data when the 
Reynolds number was greater than design. The model did 
not account for relaininarization of a turbulent boundary 
layer. The pressure gradient parameter, Ji, was calcu- 
lated. For the rear half of the pressure surface at the de- 
sign Reynolds number and below, the I\ value was greater 
than 3 x 10 -6 everywhere. At these two lower Reynolds 
numbers, the data did not indicate transition. Reference 
calculations were run for purely laminar flow at the design 
flow conditions. The results on the pressure surface indi- 
cate that a minimum Nusselt number of 250 was calculated 
at s & —0.15, similar to the calculated minimum shown 
in Fig. 6. The purely laminar predicted heat transfer then 
increased almost linearly to Nu — 600 at the pressure sur- 
face trailing edge. The experimental heat transfer is higher 
than was predicted for the purely laminar boundary layer. 
In effect it is a buffeted laminar boundary layer. 

The predictions in Fig. 6 show the consequences of ne- 
glecting the effect of turbulence augmentation in the lam- 
inar region. While the heat transfer in the leading edge 
region is now well predicted, previous comparisons (Giel 
et ah, 1999) needed augmentation in this region to agree 
well with the data. Without the augmentation, the increase 
in predicted heat transfer caused by transition near s = 0.7 
is more clearly seen at 25%' arid 50%: span. Excellent agree- 
ment with data is seen at 15% span for s > 0.5. With the 
laminar augmentation model turned on, intermitt encv in- 
formation from the CFD code indicated a transition start 
location of s = 0.2. In the transition start model, an 18% 
decrease in Reynolds number is equivalent to a 26% decrease 
in inlet turbulence intensity. Prior to transition, the suction 
surface heat transfer is now lower than the data. The pre- 
dicted spanwise heat transfer variation in Fig. 6 is similar 
to that shown in Fig. 5. On the uncovered part of the suc- 
tion surface the flow is turbulent, and the same modeling is 
used. 
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The comparisons shown in Fig. 7 at the higher Reynolds 
number show better agreement with the data for the rear 
portion of the pressure surface than did the comparison 
shown in Fig. 5 at the design Reynolds number. The over- 
prediction due to the augmentation model is the same for 
both figures. Again, the predicted midspan heat transfer on 
the uncovered portion of the* suction surface is lower than 
the data, and the 15% and 25%: predictions are higher than 
measured. The spanwise average predicted heat transfer 
is higher than measured. The experimental data in Fig. 7 
show an increase at s « 0.3 that is consistent with ear- 
lier suction surface transition. The data of Fig. 5 indicates 
that for the design Reynolds number transition occurred 
much later, at s « 0.7. For Case 3, intermittency infor- 
mation from the CFD code indicated that the predicted 
suction surface transition start occurred at about, s — 0.2, 
approximately where early transition effects are evident in 
the experimental data. 

Even though it is not shown in Fig. 8, the pre-transition 
heat transfer with augmentation over-predicted the data by 
the same ratio as shown in Figs. 5 and 7. The increase in 
Reynolds number did not affect the relative agreement. Fig- 
ure 8 shows that without the augmentation, the agreement 
is good in the leading edge region. At this high Reynolds 
number the pressure surface heat transfer is well predicted. 
The predicted suction surface transition is much closer to 
the leading edge than is seen in the data. Transition was 
predicted to begin in a favorable pressure gradient region. 
In favorable pressure gradients the location of the start of 
transition is very sensitive to the inlet turbulence level. The 
momentum thickness Reynolds number grows slowly with 
distance, and the freest ream turbulence decreases with ac- 
celerating flows. While the data gives evidence of transition 
beginning at s near 0.2, the heat transfer increase is much 
less than the predicted increase. 

Comparing the predictions in Figs. 9 arid 10 with those 
in Fig. 5 show almost identical results. Even though the 
pressure ratio changed, the exit isentropic Reynolds number 
did not. The inlet Reynolds number changed only slightly. 
Consequently, there was no mechanism to alter the pre- 
dictions. The data are almost identical, and the ratio of 
predicted to measured heat is the same as in Fig. 5. 

The predictions in Figs. 11 and 12 are with no augmen- 
tation for the laminar region. The leading edge heat trans- 
fer is well predicted. The midspan data show that, even 


though the inlet Reynolds number was only 11% greater for 
the positive incidence case in Fig. 12, the start of suction 
surface transition moved significantly closer to the leading 
edge. Comparing midspan measurements in Figs. 11 and 12 
show a small, but noticeable, heat transfer difference just 
downstream of the leading edge on the pressure surface. 
The variation is due to the different incidence angles. The 
predictions show a similar variation in the minimum values 
near s = -0.13. 


SUMMARY AND CONCLUSIONS 

Detailed aerodynamic and heat transfer measurements 
and predictions were given for a power generation turbine 
rotor under engine specific conditions. The effects of varia- 
tions in Reynolds number, exit pressure ratio, and incidence 
angle were quantified. The primary effect of Reynolds num- 
ber variations in the range of —20%. to +40% was to move 
the locations of laminar-to- turbulent transition. Variations 
in isentropic exit pressure ratio over a range of —10% to 
+ 10% were found to have little effect on the blade heat 
transfer. Positive incidence of 2 deg had almost no effect 
on blade heat transfer while 2 deg of negative incidence re- 
sulted in what could have been the start of a small pressure 
surface separat ion bubble just downst ream of the stagnation 
region. In the fully laminar or fully turbulent flow regions, 
the data agreed well with the appropriate scaling laws. 

The good spatial resolution due to the large scale and 
the liquid crystal measurement technique allowed the sec- 
ondary flow effects to be clearly quantified. These secondary 
flows significantly increased suction surface heat transfer 
rates near the endwalls. The data also quantified the three- 
dimensional impact of laminar- to-turbulent. transition on 
blade heat transfer. Excellent midspan symmetry was ob- 
served for all cases. 

Comparing computational and experimental results il- 
lustrated regions of good agreement and regions where mod- 
eling improvements are needed. This data set showed, for 
example, that improvements to the laminar augmentation 
model are needed, even though this model worked success- 
fully for previous studies from the same facility. 
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